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I The quasi-equilibrium or maximum entropy approximation is applied in order to derive constitutive equations 

. from kinetic models of polymer dynamics. It is shown in general and illustrated for an example how canonical 

distribution functions are obtained from the maximum entropy principle, how macroscopic and constitutive 
equations are derived therefrom and how these constitutive equations can be implemented numerically. In 
addition, a measure for the accuracy of the quasi-equilibrium approximation is proposed that can be evaluated 
while integrating the constitutive equations. In the example considered, it is confirmed that the accuracy of the 
OO ' approximation is increased by including more macroscopic variables. In steady elongational flow, it is found that 

04 , more macroscopic variables need to be included above the coil-stretch transition to achieve the same accuracy 

as below. 
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C3 . I. INTRODUCTION 

4— > 

In this paper, the quasi-equilibrium or maximum entropy approximation is applied in order to derive constitutive equations 
from kinetic models of polymer dynamics. It is shown in general and illustrated for an example how canonical distribution 
S functions are obtained from the maximum entropy principle, how macroscopic and constitutive equations are derived therefrom 
and how these constitutive equations can be implemented numerically. In addition, a measure for the accuracy of the quasi- 
^ equilibrium approximation is proposed that can be evaluated while integrating the constitutive equations. 
Q Polymer dynamics is a fascinating subject with important applications to polymer processing and turbulent drag reduction. 

O Very recently, direct observation of the dynamics of individual DNA molecules in shear and elongational flow have become 
possible ^ Ip. The experimental results provide useful information to improve our understanding of polymer dynamics in 
' a flow. Since deformed polymer molecules cause stresses that are responsible for the macroscopic non-Newtonian rheological 
^ behavior of polymer solutions, polymer dynamics in the presence of a flow field is also of great interest from an application point 
■rj" of view. 

OO I Simple mechanical models are very helpful in this respect. Comparison of simulation results with experiments have shown 
^n that kinetic models of polymer molecules describe very well the rheological behavior of polymer fluids U as well as the 
polymer conformation [||| in homogeneous shear and elongational flow. However, the use of these kinetic models is too time- 
consuming to be used in numerical flow calculations of polymer solutions. Therefore, simplified constitutive equations have 
been proposed in the literature that describe stresses that the polymers exert on the flow for a given deformation history (see 
"j--; for example |Q| and references therein). So far, however, no satisfactory constitutive equations are available that give accurate 

predictions under various flow conditions. 
^ Kinetic models of polymer dynamics may serve as starting point for the derivation of constitutive equations. In general, 

I ' constitutive equations cannot be derived exactly from kinetic theory due to the so-called closure problem. In addition, closure 
■ O approximations have serious effects on the behavior of the model, and sometimes even lead to unphysical predictions 

Recently, the idea of using a set of canonical distribution functions (CDF) for obtaining closures to the kinetic equation have 
been proposed in Motivated by the ease of computation, very special choices for the CDF are considered in p^ . 

In the present work, CDF are introduced using the maximum entropy principle. This approach, also known as quasi- 



equilibrium approximation, is also used successfully in other branches of statistical physics [ |1 1| , 12 1. The resulting, quasi- 
k> equilibrium approximation to the dynamics, by construction, shows several desirable features like conservation of the dissipative 
nature of the dynamics and conservation of positive-definiteness of distribution functions. We show how this approximation can 
be used to obtain closed constitutive equations for polymer solutions. 

We also present an algorithm for the numerical implementation of the constitutive equations. In addition, a measure for the 
accuracy of the approximation is suggested that does not require solutions to the kinetic model and can therefore be used while 



polymers. See also 1 14 1 for more details 



solving the constitutive equation. In the companion paper [|13|], this approach is applied to the kinetic theory of liquid-crystalline 



■ ilg@physik.tu-berlin.de 



2 



This paper is organized as follows. Kinetic models of dilute polymer solutions are reviewed in Sec. II. The maximum entropy 
principle is employed in Sec. |l|in order to obtain canonical distribution functions (CDF) for polymer solutions. Some properties 
of the CDF are also discussed. Closed form constitutive equations are derived using CDF with special emphasize on so-called 



conformation tensor models. A numerical integration scheme for the constitutive equations is presented in Sec. IV 



measure for the accuracy of the approximate description of polymer dynamics with CDF is proposed. In Sec. VI 



In Sec. a 
the derivation 



of constitutive equations using CDF together with the check of the accuracy is illustrated in detail for a kinetic model of polymer 
dynamics. Finally, some conclusions are offered in Sec. |VII[ 



II. KINETIC MODELS OF DILUTE POLYMER SOLUTIONS 



Let a linear, flexible polymer be represented by N beads. The beads themselves represent many monomers, and are large 
compared to the solvent molecules. The distance between the beads is considered to be large enough so that it allows the 
random walk statistics for the neglected segments in between, but close enough to maintain these segments in equilibrium. The 
connectivity of the chain therefore is guaranteed by an entropic potential associated with the random walk statistics between the 
beads. Ideal random walk statistics corresponds to a harmonic potential, giving raise to linear forces acting on the beads. The 
forces acting on the beads are often pictured as springs, and the kinetic models are termed "bead-spring models". 



A. Kinetic Equations of Polymer Dynamics 

Equations of motion for bead-spring models can be found in most of the textbooks on polymer kinetic theory (see e.g. 
JT^ , [1^). Let ri,i — 1, . . . , N, denote the position vector of the i-th bead. The internal conformation of the chain is described 
by the connector vectors Qi = r^+i — r^, i = 1, . . . , — 1. In the classical approach, the polymer chain is exposed to a 
given homogeneous flow field v(r; t) = vo(t) + K{t) ■ r, where K{t) = (Vrv(r; t))^ denotes the velocity gradient. Usually, 
incompressible flows, Tr k = 0, are considered. The configurational distribution function at time t is assumed to depend on 
the internal conformation but to be independent of the location of the center of mass, — 'I'({Qi}; t), where {Q, } denotes the 
set of all connector vectors {Qi} — {Qi, . . . , Qat-i}. The famous diffusion equation for 5* reads 

-^^JW (1) 

with 

N-l N-1 

i=l = l 

where Vq is the derivative with respect to Q. The last term describes diffusion due to the collisions of the solvent molecules 
with the beads. Collisions keep each bead in thermal equilibrium with the solvent. Temperature, Boltzmann's constant and 
friction coefficient of the beads are denoted by T, fee and (, respectively. The function U ({Qi}) is the dimensionless potential 
of forces on the beads so that Fi({Qi}) = k-QTVc^JJ is the force acting on the i-the bead. It is generally assumed that only 
adjacent beads interact, the interaction being isotropic, so that U can be written as 

N-l 
i=l 

Long range interactions along the chain, such as those leading to excluded-volume effects might require some extra treatment. 
The diffusion matrix D^j can be represented as = D{Af^l+C,^lij ), where the diffusion constant D is given hy D — 2fcBr/C- 
The unit matrix is denoted by 1, and the (A^ — 1) x [N — 1) Rouse matrix is 

AR:=2%-(4,+i+,5,,j_i). (4) 

The tensors Clij ({Qi}) describe hydrodynamic interactions of bead i and j as the movement of the beads disturbs the flow field. 
Typically, the hydrodynamic interactions are described by (regularized) Oseen-Burger tensors [|l6|]. The first term on the right 
hand side of Eq. (^ accounts for the hydrodynamic drag the beads experience while moving through the flow. It is assumed that 
the drag force is proportional to the difference of bead and solvent velocity. 

The number density of polymers rip — (1) is conserved by the dynamics (^. Here and in the sequel, brackets (•) denote 
averages with respect to ^, (•) = ^dfiQ • *({Q.}) with d = i{N - 1). Eqs. (|l])) and (||) describe an ensemble of non- 
interacting bead-spring chains in the presence of a flow field. Therefore, these models describe the dynamics of dilute polymer 
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solutions. The Rouse model, for example, is obtained when hydrodynamic interactions are neglected, f2y = 0, and harmonic 
potentials V^(Q^) = HQ,^/ {2k'QT) are considered, H denotes the spring constant. However, harmonic potentials are unrealistic, 
since they allow infinite stretching of the chain in strong flows. Warner's finite extensible nonlinear elastic potential (FENE), 

y(Q2) = _(6/2)hi(l-QVQLx), (5) 

exhibits a logarithmic divergence at the maximum length Q^iax, and therefore prohibits further stretching Jl7|]. The finite 
extensibility parameter, h — HQl^^^/ (kBT), is a dimensionless measure of the maximum extension. The Hookean spring force 
is recovered in the limit 6 — > oo. 

We restrict ourselves here to finite extensibility effects with or without hydrodynamic interactions. Other effects such as 
excluded volume, internal viscosity, friction coefficient varying with extension or anisotropic frictional drag can also be included 
in bead-spring models (see [|l5]| and references therein). 



B. Entropy Production 

In the absence of flow, the distribution function relaxes to its unique equilibrium state, = jZ, due to the dynamics 
(|l|) and (^, where Z is the normalization constant. The Boltzmann entropy functional, 

= -fcB J d^Q *({Q.}) ln[vI/({Q,})/vI/^<i({Q,})], (6) 

is strictly non-decreasing during the relaxation due to the dynamics (|l]) and for /« = 0. The entropy production S is defined 
as the time evolution of the entropy functional (0), 



where 5/ 5"^ is the Volterra functional derivative. According to Eqs. and (Q), the entropy production can be split in two parts, 
S ~ S^^ + S^, corresponding to the first and second term on the right hand side of Eq. (^, respectively. It is easily verified that 
S is non-decreasing in the absence of flow. 



c-d _ Kb 



if the diffusion matrix D is positive semi definite. In the presence of a flow field, the entropy exchange 5'' of the polymer and 
the solvent subsystem is given by 

5" = -i/^:xP (9) 

which is not sign-definite. In Eq. (^, the polymer contribution to the stress tensor tP has been introduced which is given by 
Kramer's expression [ p3| , 

TP[m]^[N -l)npkjiTl-kBTY, J d''Q■^{{Cl^})^^VQ,U{{Q,}). (10) 



III. CANONICAL DISTRIBUTION FUNCTIONS IN POLYMER DYNAMICS 



Within kinetic theory, the state of the polymer at time t is described by the configurational distribution function ^({Q;}; t) 
and Eqs. (|l]), (^ and ( p^ specify the dynamic and viscoelastic behavior of dilute polymer solutions for a given homogeneous 
flow field. In viscoelastic flow calculations, the velocity field has to be found numerically together with the dynamics of the 
polymers. The combined simulation of Navier-Stokes equation with the kinetic equation ([l|) is however very time-consuming. 
Therefore, there is an enormous interest in deriving closed constitutive equations from the kinetic equation. In general, this 
derivation cannot be done rigorously due to the so-called closure problem and certain approximations have to be invoked. 
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A. The Closure Problem 



It is assumed that the macroscopic behavior of the polymers can be described by a set of n macroscopic variables, which are 
defined as functionals of the distribution function, Mk = nfe[^'], k ~ 1, . . . ,n. We consider macroscopic variables that are 
linear functionals (moments) of the distribution function "if 

Mk{t)^ Jd''Qmki{Q^})'i'i{Q^}■,t). (11) 

The set of macroscopic variables need not be identical with the components of the stress tensor tP. We only require that the 
stress tensor (jl^) is determined by the values of the macroscopic variables M — {Mi, . . . , Af„}. From the chain rule of 
differentiation and Eq. (|]), the time evolution of Mk is given by 

Mk = J d'^Q mfe({QJ) J(^') (12) 

which do not constitute, in general, a closed set of equations for the macroscopic variables. In what will follow, we present a gen- 



eral approach to obtain closures of Eqs. ( 12 1 based on the maximum entropy principle and corresponding canonical distribution 
functions. 



B. Maximum Entropy Principle and Canonical Distribution Functions 

The reduced or macroscopic description of the system is reasonable if a separation of time scales of the relaxation times 
holds. The equilibrium distribution function ^''^^ = e^^ /Z is obtained from the maximum of the entropy functional (^ subject 
to fixed normalization, 1 = j d'^Q ^'({Qi}). The so-called generalized canonical distribution functions 4** are obtained by the 
maximum entropy principle under fixed constraints [ |l2[ [l8| , [l9[ ^0[ |, 

^ max, J d^'Q mfc({QJ)*({Q,}). (13) 

Solution to the problem Eq. ([l^), if it exists, is unique for concave entropy functionals like (^ and linear constraints given by 
Eq. (11). For the present case the solution to Eq. ( pj[ ) is given explicitly as 

n 

M'*({QJ) = *-'i({Q,}) exp[^ AimiiiQ,}) + Aq] (14) 

1=1 

and is called canonical distribution function (CDF) (or generalized canonical state or quasi-equilibrium distribution function or 
maximum entropy state). The Lagrange multipliers, A = {Ao, Ai, . . . , A„}, are determined from the constraints 

„ n 

7ip = yd'^Q*°'i({Q,})exp[^A,mi({Q,})+Ao], (15) 



and 

d'^Q mk{{Q^}W{{Q^}) cxp[^ A,mi{{Q,}) + Aq]. (16) 

;=i 

In most cases, the dependence A(M) or M(A) is not known explicitly. The manifold ft* that is spanned by when the 
macroscopic variables take on all admissible values is known as quasi-equilibrium manifold or canonical distribution function 
manifold [mH]. 

The macroscopic entropy 5** (M) is defined as the entropy given by Eq. (||) evaluated on the CDF (|l^, S* (M) = S[^*]. 
The macroscopic entropy S* allows to interpret the Lagrange multipliers as the conjugate or dual to the macroscopic variables. 



(17) 
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C. Macroscopic dynamics 

In what will follow, the canonical distribution function manifold H.* is used to obtain the macroscopic description of the 
polymer dynamics. We assume that arbitrary initial states relax quickly towards the manifold fl* and further evolve in the 
neighborhood of In this case, the CDFs that span the manifold can be considered as representative states in the sense that 
moments of the solution to the kinetic equation (|l]) are approximated accurately by the corresponding moments evaluated with 
the CDFs. The validity of this assumption is closely related to the separation of time scales mentioned above and the choice of 



the macroscopic variables. This assumption can be interpreted within the context of the method of invariant manifold |21 1 as 
O* being a meaningful approximation to the relevant invariant manifold. A measure of the accuracy of this approximation is 
presented in Sec. ^ For improvements on this approximation see [ pi] , p2| ]. We do not provide any rigorous proofs for the results 
obtained in the sequel. However, the ability of measuring the accuracy of the approximation is a constructive method to gain 
confidence in the approximation scheme proposed here. 

Having identified the manifold fl*, the macroscopic time evolution is defined as 

Mfe* = J d^Q mfc({Q,}) J(**(M, {Q,})), (18) 

which is the dynamics of the macroscopic variables ( [T^ evaluated on fi* . Inserting the right hand side of Eq. (||) for J, Eq. ( [l8| ) 

becomes 

N-l n 



1=1 1=1 

where (•)^ denotes averages with respect to the CDFs 5**, Eq. (p^), and AI is a positive semi-definite matrix, 

1 f 

Mki = ^J2 J ^""Q (VQ.™fe) • • (Vq^mz)- (20) 

i J = l 

Eq. ( [l9| ) together with Eqs. (|l5|), dl^ ) and ( po| ) represent the set of closed equations describing the macroscopic dynamics. Since 
the CDFs are assumed to be representative distributions, the constitutive relation is obtained by evaluating Eq. ( p^ on D,*, 
T-P*(M) = -rP[^'*]. In the present case it reads 

N-l n 

tP*(M) = -ksT J2 5I(Q'^Q."^fe)AAfc- (21) 

i=l fe=l 

From Eqs. (jl^) and ( p^ it is found that the macroscopic entropy production S* is given by 

1 " 

S* = --k: xP* + ^kMkiAi- (22) 

k,l = l 

and could have been obtained by evaluating Eqs. (||) and (^ on the CDF (|l^. In the absence of flow, the macroscopic entropy 
production (|2|) is positive semi-definite as is the underlying kinetic model, Eq. (8^ This illustrates a general property of the 
quasi-equilibrium approximation that it conserves the type of dynamics (see, e.g. 122]). 



D. Conformation Tensor Models 



In the present section, we specialize to conformation tensors as macroscopic variables. The corresponding conformation tensor 
models are frequently used in the literature p^]. The conformation tensors cr are the covariance matrices of the distribution 
function 

Ty = (Q,Q,), = ,N-1. (23) 

In our notation, conformation tensor models are defined by the special choice Mk cTij and rrik QiQj in Eq. (pi]). 
Corresponding Lagrange multipliers are denoted by Ay . 

For conformation tensor models the time evolution equation ( p^ ) reads 

n 

a-ij = K ■ (Tij + <Tij ■ k'^ - Y Mi-jki ■■ ^ki, (24) 

k,l=l 
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where the matrix ( |20| ) now reads 

M,jki = ^ (Q^DJfeQ^ + Q,Dj7Qfe + QjD,:fcQi + Qjf) aQk) ^ . (25) 

In Eq. ( p4] ) one recognizes the upper-convected derivative, while the last term describes relaxation mechanisms. 
For the present choice of macroscopic variables, the stress tensor ( pTj ) takes the form 

N-l 

rP = -2fcBT ^ cr„ • A,,. (26) 



IV. NUMERICAL INTEGRATION SCHEME 

In the maximum entropy principle (pj[), Lagrange multipliers A occur naturally to satisfy the constraints. In principle, the 
Lagrange multipliers can be expressed as functions of the macroscopic variables using the constraints (|5]) and (|l6|). In general, 
however, no explicit expressions A(M) can be found. For these cases, we here propose a method that does not require knowledge 
of the functions A(M) 



A. Dual Variables and Their Dynamics 

A change of the independent variables from the set M to A is done by Legendre transform. The Legendre transform of the 
macroscopic entropy S* (M) is the generating function G(A) = S* (M) + J2k ^kM^, which can be written as 

/n 
d^Q *^<i({Q,})exp[^ A;mz({Qa)]. (27) 

The macroscopic variables and higher order moments can be obtained by derivatives of G with respect to A, 

Mfc(A) EE (mfc)^ - (28) 

and 

CkiiA) = {mumi)^ - (mfc)^ {mi) ^ = (29) 

The Lagrange multipliers A are now considered as primary independent variables that parameterize the CDFs (|l^ and that 
determine the values of the macroscopic variables (28 ). From Eqs. (^Sj) and (^9|) it is found that the time evolution of A is given 
by 

n 

^k^Y.'^ki\A)Mi{A), (30) 

1=1 

where Affc(A) is given by Eq. (|l8|). Eq. (|o|) is the closed macroscopic equation written in dual variables. Note, that the 
constitutive relation (|T]) is already given as a function of A. 



B. Numerical Integration Scheme 



The reformulation of the macroscopic dynamics ( |19[ ) in terms of the dual variables A described in Sec. IV A is suitable for 
numerical implementation. The Lagrange multipliersA now play the role of independent dynamic variables, instead of M_. To 
advance given values A{t) at time t to their values A(t + r), with small time step r, the following first order integration scheme 
is proposed |24|: 



1 . The new values of the macroscopic variables M(i + r) are found from the values A{t) by 



Mkit + t) - Mkit) 



d'Qm,({Qa)J(**(A(t),{Q,;})), 



(31) 



where J(** (A(i), {QJ}) denotes the right hand side of Eq. (^, evaluated with the CDF ([l4|) with A{t) the actual values 
of the Lagrange multipliers. 
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2. The matrix Cki, defined in Eq. (|29|), is evaluated from 

C'ki{t) = (TOfem,)A (t) ~ Mkit)Mi{t), (32) 
where {mkrni) ^ is calculated from Eq. (pTj), where is now replaced by nikmi. 

3. The n X n matrix C(<) is inverted numerically to give C~^{t). 

4. The values of the Lagrange multipliers, A{t + r), are given by 

Akit + T)-Ak{t)_j^^_,^^^Mk{t + T)-Mk{t) 
1=1 



This concludes one time step of integration. The integration scheme has to be supplemented by initial conditions A(0). The 
special case of equilibrium initial conditions corresponds to A(0) = 0. 

Note, that evaluating the matrix ( p9[ ) requires moments of which are of higher order than the macroscopic variables them- 
selves. If n denotes the number of macroscopic variables, the numerical integration scheme requires in total n(n + 3)/2 integrals 
per time step to evaluate M and C. Due to the symmetry of C, this number is of order n^/2 rather than n'^. The evaluation 
of Yfe may require additional integrals. The evaluation of all these integrals with standard numerical methods might be time- 
consuming, especially for high-dimensional integrals. It is demonstrated in [ ^5| ] that under certain circumstances these integrals 
can be evaluated efficiently by adapting methods of fast Fourier transformations. 

V. ACCURACY OF APPROXIMATION BY CANONICAL DISTRIBUTION FUNCTIONS 

The reduced description outlined in Sec. |l| obviously imposes restrictions or approximations to the underlying kinetic model. 
Therefore, the accuracy of the reduced description needs to be discussed. The convergence of arbitrary initial conditions to 
the manifold fl* is of obvious importance, but is in general difficult to study. For a concrete example, the convergence or 
non-convergence is investigated in [p6|]. 

CDFs are in general not solutions to the kinetic equation i.e. the canonical distribution function subspace fl* is not 
dynamically invariant. The dynamic variance A is defined as 

A(M) = J(vI/*)-f]|^M,, (34) 

k=l ^ 

which is the difference between the microscopic time evolution evaluated on Vl* and the time evolution of the CDFs due to 
the macroscopic dynamics. The dynamic variance A can be considered as 'defect' of the macroscopic description. While the 
method of invariant manifold [|T]] uses A to obtain improved manifolds, we will later use | A| as a measure of the accuracy of 
the approximation by CDFs. The norm | • | used later is related to A^, which is the difference of the dynamics of moment 
due to the microscopic and macroscopic dynamics, 

/d 
(fQ meA = - [(m^mj)A - {mp) ^{mj) A\{C-^)jkMk. (35) 

A^ vanishes identically for £ < n as it has to. For £ > n, the quantities {mg)/^ are not included in the list of macroscopic 
variables ([l3|). If these quantities are defined such that the set {M, {mi)^} is linear independent of M_, Eq. ( ^ provides 
valuable information on the dynamics of variables that are not included in the list of macroscopic variables. 

Inserting in Eq. (Q) the right hand side of Eq. for J and Eq. (|l^ for ^!* together with the time evolution of the macroscopic 
variables (|9]) one obtains Am = Am with 

7V-1 N-1 



AM = -TrK-^Q, •/«t.(-Vq.[/ + ^°)- ([VQ."^fe]Q^)ACfe (36) 

i—1 i,k—l 
N-1 N-1 



2 

i.j — 1 i.j,k=l 



In Eq. (^6|), we have defined = ^^j, A^Vq^mfc and — din ^'*/aAffe. For the special case of conformation tensors as the 
macroscopic variables the expression ( ^6| ) simplifies on the replacement M QQ. If in addition one specializes to certain 
mean-field potentials, then the expression for the variance becomes most transparent. This is demonstrated in the sequel. 
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A. Gaussian Manifolds of Conformation Tensor Models 

In this section we consider general potentials of the form 

N-l 

U{W^,})=J2Vi^^'^^^) (37) 



for the kinetic equation (|l]). The restriction to these potentials is very important. They can either be viewed as mean-field 
approximations to an underlying potential of the form (||) or as models in their own right. Concerning the first interpretation, a 
possible derivation together with systematic corrections to this approximation is illustrated in 

For conformation tensor models, the equilibrium states are isotropic Gaussian distributions if the potentials are of the form 
(|3^). An important property of conformation tensor models for these potentials is that also the quasi-equilibrium states ( p^ are 
Gaussian. In this particular case, Eq. ( ^3| ) can be evaluated by Gaussian integration and subsequently be inverted to give the 
Lagrange multipliers as explicit functions of the macroscopic variables 

Ay({<T.,}) = y'(Tr^,,)l% - li^^-'h- (38) 

Inserting the functional form of the Lagrange multipliers (38 1 and the potential ( p7| ) into the quasi-equilibrium distribution 
function we obtain the standard representation of the symmetric Gaussian distribution functions 

JV-l 



V(27r)'^ deter ^ 2 /j^^ 
Up to an irrelevant constant, also the quasi-equilibrium entropy S*{(t) is found from Eq. (|3^) as explicit function of cr. 



(39) 



S*iW^J}) = - In dot a- - C/({^,,}). (40) 



The first term on the right hand side of Eq. ( |4-0| ) originates from the normalization of ^cr while the second term contains the 
entropic spring potential. 

B. Variance of Conformation Tensor Models 

For potentials of the form (p7|), the time evolution equation for the conformation tensor cr is obtained from Eqs. (p4|), and 
( |3^ ) and reads 

JV-l 

O-y = ft • CTij + (Tij ■k'^ V'(Tr(Tkk){Q,t>jk ■ Qfe + D.fc • QfeQj)A + (Dy)A, (41) 

fc=l 

where use has been made of the decomposition of moments for the Gaussian distribution functions (see e.g. Appendix A of 
p^]). It is not surprising, that Eq. (^Ij) coincides with the result of the Gaussian approximation used in p7[|, the averages 
being performed with the Gaussian distribution (|3^). Eq. ( |4l| ) is an obvious generalization of Eq. (5) in [|7|| for the anharmonic 
potentials of the form (|3^). In Eq. ( |4l| ) the "incompressibility" condition V • D = for hydrodynamic interaction tensors is 
assumed. 

The variance of the QEA (|6|) can be calculated for the kinetic equation (jl]) by Gaussian integrals. The result is 

JV-l N-l 

Ao- = ^o-h 51 Tr(T„ -(o--i),,)+ ii'^^'U-Qkf ■T,,-{{a-^)ji-qi)l (42) 

where we have introduced the tensor Ty — + Tjj, with the contribution of the Brownian part as T^^ — T)ij — (D^ }a and 
the spring contribution as T^^ , 

AT-l 

T'j = Y V'{Tl■CTkk)[{^D^l ■ QiQj)a - t)^l ■ CTij]. (43) 
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We will see in a moment why the flow term does not contribute to Ao-- 

In the absence of hydrodynamic interactions or for pre- or self-consistently averaged hydrodynamic interactions, the diffusion 



matrix D is independent of {Qi}. Upon inspection of the variance (42), one notices it immediately that in these cases the 
variance Ao- vanishes identically. In these case, the drift term in the diffusion equation (Q) is linear in Q and the diffusion term 
is independent of Q. In this case, Gaussian distributions provide exact solutions to the diffusion equation [|2^. Examples are 
the Rouse and Zimm models but also models that incorporate finite extensibility effects in the Peterlin approximation. From the 
above discussion we learn that in this case the Gaussian distributions constitute the exact invariant manifold of the dynamics (|l]) 
and Eq. (^l]) represents the time evolution equation on the invariant manifold. Since there exist many solutions to Fokker-Planck 
equation the relevance of the Gaussian solutions is discussed in the next section. 

For general hydrodynamic interactions the variance does not vanish and the quasi-equilibrium approximation (|9|) is no 
longer exact. It may serve as initial approximation as, for example in [^. Possible improvements on the QEA can be found in 

VI. ILLUSTRATION 

As we have already mentioned, the choice of macroscopic variables is not unique. Very different choices of macroscopic 
variables are possible, and correspond to different sets of microscopic states accessible to the macroscopic description. In the 



literature, conformation tensors are frequently used as macroscopic variables. In Sec. V B it is shown for a particular class of 
models, that the quasi-equilibrium manifold for these variables is the invariant manifold. In case the invariant manifold is stable, 
one observes a fast relaxation to the invariant manifold followed by a slow evolution on the manifold. 

In general, the quasi-equilibrium manifold is not invariant and improvements on this approximation become an important 
issue. If the set of macroscopic variables is maintained, successive improvements on the QEA can be found along the lines of 
the method of invariant manifolds [plt|. Another option is to extend the set of macroscopic variables and apply the QEA for the 
larger set of variables. Here, we illustrate the latter option for the kinetic model of polymer dynamics considered in [p9||. 

In the context of polymer dynamics, the problem of the choice of macroscopic variables is encountered frequently. Recent 
experiments of Chu and co-workers ^ clearly demonstrate the dependence of the polymer dynamics in flow on the initial state. 
A possible classification of initial states into dumbbell, half-dumbbell, folded, coiled and kinked conformations is proposed in 
1^]. The failure of the macroscopic description in terms of a single conformation tensor has attracted considerable interest and 
led de Gennes to the notion of 'molecular individualism' [^0|]. 

In order to deal with this situation, the macroscopic description requires an extended set of variables. In |l^], the addition 
of one extra variable associated to the fourth moment of the distribution function is proposed. In [^, [l^], also a closure approxi- 
mation for this set of variables is proposed. In this closure, however, a penalty function has to be introduced ad hoc in order to 
obtain physically admissible results. 

Here, we study the macroscopic description for an extended set of variables in the QEA. As demonstrated in [p^, the QEA 
is well-behaved for any choice of macroscopic variables. In the limit where the number of macroscopic variables goes to infin- 
ity, the microscopic description is recovered since no variation over microscopic states is performed in the maximum entropy 
principle and therefore no reduction of description is associated with the QEA. We illustrate here for a finitely extensible dumb- 
bell model that the macroscopic description converges fast with increasing number of macroscopic variables. We consider the 
one-dimensional model proposed in [|9[, with the spring force 

FiQ) - HQi tan(Q/gi), for - Q,nax <Q< Qmax (44) 

where Qi — 2(5niax/7i' is determined by the maximum extension Qmax and H denotes the "spring constant" in the sense 
F{Q) HQ for Q — > 0. The force ( pj ) approximates the inverse Langevin function very well and allows the exact calculation 
of the spectrum of relaxation times of the associate Fokker-Planck operator The result is 

A„ = 2Ah [?i + nV6i]~\ n = l,2,... (45) 

where Ah = C/ i'^H) is the time constant of the Hookean dumbbell and bi = HQ\/k-QT a dimensionless finitely extensibility 
parameter [E9||. Also the exact eigenfunction have been found in [^91] to be given by — ^n^V<i with 



= C^ala\ . . . al_, cos^+'^/^Q/Q^) (46) 

where aj^ = — Vq + /q,^ are generalized creation operators. C„ denote normalization constants and fa^. = atan{Q/Qi). 

The first eigenfunctions read: *o = *cq, *i = ^^i sm{Q/Qi)'^cq, *2 = A^2[(fei + l)sm'^{Q/Qi) - coa'^ {Q / Q i)]'^ cq with 
Ni = bi + 2 and A^l = (6i + 4)/[2(6i + 1)]. Note, that the eigenfunctions are orthonormal, ('!'„, ^'m)s = Snm, with respect 
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to the scalar product 

(5,/»)s= / dQ^-^{Q)g{Q)h{Q). (47) 

"'-Qmax 

As macroscopic variables, choose the first n even dimensionless moments of the distribution function = (Q^'^)a^ k — 
1, . . . ,n and Q — Q/Qi is the reduced extension. The benefit of considering a one-dimensional model is of course that there 
is no tensorial structure associated with these moments. The QEA for this choice of macroscopic variables is of the form ( |l4| ) 
and reads 

n 

^X(Q) = (l/-2cq) cos^^ (Q) exp [J2 AjQ'^ + Ao]. (48) 



The dynamic equations for the macroscopic variables found from Eq. ( |19| ) are 

Mfe = 2kKMk + k{2k - l)Mfe_i/(6iAH) - kRk/Xu (49) 



where Rk — (tan((5)(5^'^ ^)a. The polymer contribution to the stress found from Kramer's expression ( pX]| ) is t'p /[upk^T) — 
l-hRi. 

For given values of the n Lagrange multipliers at time t, A{t), the corresponding values of the n macroscopic variables are 
found from 



2„Ao(t) .7^/2 _ . " . 

Mk{t) = / dQ Q2fc cos'^ (Q) exp[^ A{t)Q^'] (50) 



Thanks to the finite extensibility of the model, the integration in Eq. (^0|) is performed over a closed interval only. In addition, due 
to the symmetry of the integrand, the integration is performed over the positive range of Q only. The matrix Cki — (Q^'''^'' ) a is 
obtained by the n{n + 1) /2 corresponding integrals. In this particular case, some of these integrals are identical to components 
of M. Inversion of the matrix C is facilitated by the fact that it is a symmetric, positive semi-definite matrix. In total, one 
time step of integration requires n{n + 5)/2 one-dimensional integrals similar to Eq. (|o|). The numerical integrations were 
performed using standard integration schemes such as Gauss-Legendre or Romberg integration [pi[|. The algorithm, described 



above and in Sec. IV B was implemented and run for fixed time step r in the order of 10^ . . . 10^ Ah. 

In particular, we consider equilibrium initial conditions (Afc(i = 0) = 0) and steady flows and compare the quasi-equilibrium 
dynamics for different numbers n of ma croscopic variables. We generally expect the QEA to be a good approximation near the 



true equilibrium. As noted in Sec. [II D, the approximation ( [4-8[ ) is exact in the steady state since the one-dimensional analogue 
of the conformation tensor is kept as macroscopic variables and any flow in one dimension is a potential flow. Therefore, we 
expect (j^m to be a good approximation also near the steady state. 

In Fig~|], the polymer stress t'^ / [ripk^T) is shown as a function of dimensionless time i/An for equilibrium initial conditions 
and different strengths of the flow. The finite extensibility parameter was chosen to be 6 ~ 50. Det ails of the numerical 



integration of Eq. (^^ and evaluation of Rk and rP with the distribution function ( }48| ) are given in Sec. |IV B| . As seen from 
Fig. |l] a), we observe that for weak or moderately strong flows, k ^0.5, the analogue of the conformation tensor (n — 1) already 
provides a satisfactory description. Within the resolution of Fig. |l] a), the predictions are the same also for higher values of n. 
We observe that the description in terms of the second moment becomes worse in stronger flows, k > 0.5 Including higher order 
moments, we observe a fast convergence of the polymer stress in the macroscopic description to the result of the microscopic 
simulation. As seen in Fig. |l]b), the first three even moments already provide a satisfactory description for k = 1. 

It is interesting to trace back the difficulties of the macroscopic description in terms of second moments to the behavior of 
the probability distribution function (pdf) in these flows. Since the macroscopic description corresponds to the manifold of 
distribution functions (p8[), direct comparison between microscopic and macroscopic description is possible also on the level 
of distribution functions. In Fig. ^ the quasi-equilibrium states (^sj) for the corresponding values of the Lagrange multipliers 
are compared to the numerical solutions of the kinetic equation at time t = 2Ah and t — 8Ah- In the present case, the kinetic 
equation is not solved directly but the stochastic differential equation associated with the kinetic equation is solved numerically 
for a large ensemble of dumbbells |]l6[]. This method is well-established and is known as Brownian dynamics simulation. In this 
case, the distribution function is represented by a histogram over the lengths of the dumbbells of the ensemble. As seen from 
Fig. I a), the QEA for n = 1 gives a good approximation even for the distribution function in case the elongational flow is not 
too strong, k < 0.5. In stronger flows, the shape of the distribution function is not accurately predicted by the quasi-equilibrium 
distribution (Q) for n = 1, both, for rather early times (circles) and later times (diamonds). A much better approximation of 
the distribution function for earlier times is obtained for n = 2. Note, that for later times a characteristic peak of the distribution 
function near the maximum extension occurs indicating that most dumbbells are stretched to almost their maximum length. 
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Note, that the precise shape of that peak is not described perfectly by the quasi-equilibrium distributions even for n = 3. The 
approximation is however reasonably close to give accurate predictions for the polymer stress as seen in Fig. |l]. 

As discussed in Sec. ^ the accuracy of the approximation of the dynamics by canonical distribution functions can be measured 
by a suitable norm of the dynamic variance A. Let us reconsider the approximation (^) with n = 1. Then, the quantity A2 
defined in Eq. ( p^ ) for £ = 2 measures the difference of the microscopic and macroscopic dynamics of the fourth moment. 
For two strengths of the flow, A2 is given as a function of time in Fig. 3|. From Fig. ||one notices that A2 vanishes for early 
times i/An ^ 1 since the initial state was chosen to be the equilibrium distribution for which A2 = 0. Also for times i S> 1, 
the dynamic variance A2 vanishes, since the steady state is correctly described by the approximation (^8|). For intermediate 
times t « 1, A2 is different from zero and can be interpreted as a measure for the quality of the approximation. For k — 0.5, 
A2I remains small, less than 10^'^, while for k = 1.0, IA2I is two orders of magnitude larger This is in agreement with the 
comparison to Brownian dynamics simulations shown in Fig. ^that the approximation n = 1 performs quite well for k = 0.5, 
whereas it is much less reliable for k — 1.0. Note that | A2 1 is a measure for the difference of the microscopic and macroscopic 
time evolution. Therefore, the maximum of IA2I occurs at earlier times than the maximum difference of the moments and the 
polymer stress. 

Although the model studied in this section is only one-dimensional, we expect the basic observations to be valid also in three 
dimensions. The present approach can easily be generalized to three dimensions and it would be interesting to confirm this 
conjecture. Remember that the quasi-equilibrium distributions (Q) are exact in equilibrium and in the steady state in case of 
potential flows. Therefore, it would be particularly interesting to see whether the above results also apply in case of non-potential 
flows in three dimensions. 

It is straightforward to apply the numerical integration scheme proposed in Sec. [V B to the three-dimensional problem. 
The corresponding numerical integrations are however much more time consuming than in the one dimension. Under certain 
circumstances, these integrals may be evaluated efficiently, e.g. by fast Fourier transform algorithms [p5||. 

Based on the observations for this model, we conclude this section with the remark that the frequently used conformation 
tensor models might not perform very well in strong flows if compared to the original model. In the present example, the failure 
of the conformation tensor model is_related to the so-called 'coil-stretch transition' [Q, giving raise also to the instability of 
the invariant manifold as studied in 



VII. CONCLUSION 



In the present work, we have shown how the quasi-equilibrium or maximum entropy approximation can be applied to kinetic 
models of polymer dynamics. As a result, canonical distribution functions and constitutive equations of polymer dynamics are 
obtained. A numerical implementation of the constitutive equations is proposed together with a measure for the accuracy of 
the quasi-equilibrium approximation. The numerical implementation and the measure of the accuracy of the approximation 
is illustrated in detail for an example of polymer dynamics. It is found that the accuracy of the approximation increases with 
increasing number of macroscopic variables. In steady elongational flow, more macroscopic variables are needed in order to 
arrive at the same accuracy above the coil-stretch transition than below. 

Finally, we mention that the same methods can be applied to other kinetic models of polymer dynamics like, e.g. liquid- 
crystalline polymers [ pj[ ] . 
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FIG. 1: Performance of QEA in steady elongadonal flow with k = 0.5 a) and k = 1.0 b). Shown is the polymer stress in units of n^k^T 
as a function of dimensionless time t/Xa for different sets of macroscopic variables. Dashed: n = 1, dotted: n = 2, solid: n = 3 and 
symbols: Brownian dynamics simulations of full kinetic model. 
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FIG. 2: Performance of QEA in steady elongational flow with k = 0.5 a) and k = 1.0 b). Shown is the distribution function as a function 
of reduced extension Q/Qmax for different sets of macroscopic variables at two times. Dashed: n = 1, dotted: n = 2, solid: n = 3 and 
symbols: Brownian dynamics simulations of full kinetic model. 
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FIG. 3: Dynamic variance A2 as a function of reduced time t/An- Solid line corresponds to k = 0.5, broken line to k = 1.0. Comparison to 
Fig. [l| shows that the dynamic variance A2 can be used as a measure for the accuracy of the approximation. 



